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ABSTRACT 

The  method  of  vanishing  artificial  viscosity  is  used  to  obtain  smooth, 
large-data  travelling-wave  solutions  to  a  class  of  conservation  laws  with 
semidefinite  viscosity.  The  one^dimensional  Navier-Stokes  equations  serve  as 
an  illustrating  example.  ■ 
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§1.  INTRODUCTION 

In  this  paper  we  discuss  sufficient  conditions  for  the  existence  of 
smooth,  large-data  travelling-wave  solutions  to  certain  systems  of  nonlinear, 
semidefinite  parabolic  partial  differential  equations.  Specifically,  we 
consider  conservation  laws  which  have  the  form: 


(1.1) 


ft  *  ftr1 '  fe  (“■>  ft)  • 


where  u,  f(u)  €  and  where  A(u)  e.  RPxn  may  be  singular.  Moreover,  we 
include  in  this  study  problems  for  which  f(u)  and  A(u)  are  not  globally 
defined.  To  illustrate  our  ideas  we  apply  the  abstract  results  of  this  paper 
to  the  one-dimensional  compressible  Navier-Stokes  equations;  see  [6]  for  a 
classical  treatment  of  this  problem.  In  the  Navier-Stokes  equations  the  left 
hand  side  of  (1.1)  is  given  by: 
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where  p,  m  and  e  are  respectively  a  fluids  density,  momentum  and  total 
energy.  The  variable  P  represents  the  fluids  pressure  and  is  given  by: 

P  -  (Y  -  1)(e  -  m2/2p)  , 

for  an  ideal,  polytropic  gas.  y  >  1  is  a  thermodynamic  constant.  The 
viscosity  matrix  A(u)  is  given  by: 
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where  again  A  and  cv  are  thermodynamic  constants,  and  p  is  a  positive 
constant  often  referred  to  as  the  coefficient  of  viscosity. 


A  travelling-wave  solution  is  a  particular  solution  of  (1.1)  which  has 
the  form: 


(1.2)  u(x,t)  =  w(x  -  st)  , 

lim  w(x  -  st)  =  u^  , 
x+±“ 

for  some  scalar  constant  s.  Inserting  the  Ansatz  (1.2)  into  (1.1)  and 
integrating  we  find  that  w(t)  satisfies  the  ordinary  differential  equations 

(1.3)  -s(w  -  u_ )  +  f(w)  -  f(u_)  =  A(w)  , 

where  here  t  represents  the  travelling-wave  variable  x  -  st.  (In  later 
sections  we  use  t  to  represent  the  travelling-wave  variable.)  For  a 
trajectory  of  (1.3)  to  satisfy  (1.2)  it  is  easily  seen  that  given  u_  we  must 
have  that  u+  and  s  satisfy  the  Rankine-Hugoniot  condition: 

(1.4)  -s(u+  -  u_)  +  f(u+)  -  f(u_)  *  0  . 

(From  time  to  time  we  shall  write  s  as  s(u_,u+)  to  indicate  its  explicit 

dependence  on  the  states  u_  and  u+.)  Given  that  u_,  u+  and  s  satisfy 

the  Rankine-Hugoniot  condition,  it  is  a  natural  question  to  ask  whether  there 

is  a  smooth  function  w(t)  that  satisfies  (1.2),  (1.3).  To  answer  this 

question  when  the  viscosity  matrix  A(u)  is  singular  our  approach  is  to 

perturb  A(u)  to  a  nonsingular  matrix  AC(u)  (with  A°(u)  *  A(u))  and 

establish  the  existence  of  a  family  {wG}  of  solutions  to  the  modified 

e 

problems.  The  behavior  of  w  as  e  tends  to  zero  is  then  investigated. 

For  the  Navier-Stokes  equations  the  modified  problems  are  obtained  by 
introducing  the  artificial  viscosity 


(1.5) 
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(with  a(u)  >  0  to  be  determined),  into  the  conservation  of  mass  equation 


This  paper  is  divided  into  five  sections.  In  Section  2  we  list  a  number 


of  preliminary  assumptions.  Section  3  is  devoted  to  the  existence  of 
travelling-wave  solutions  to  (1.1)  when  the  viscosity  matrix  A(u)  is 
nonsingular,  Fingular  viscosity  problems  are  treated  in  Section  4.  Finally, 
in  Section  5  we  apply  our  results  to  the  compressible  Navier-Stokes  equations 

Acknowledgement .  We  wish  to  thank  Professor  J.  V.  Ralston  of  UCLA  in 
particular  for  suggesting  the  proof  of  Theorem  1  as  well  as  for  his  other 
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§2.  SOME  PRELIMINARY  ASSUMPTIONS 

The  global  structure  of  the  Hugoniot  locus  (that  is  the  set  of  states 
u+  that  satisfy  (1.4)  for  some  fixed  u_  and  variable  s  e  R)  has 
fundamental  importance  to  our  analysis.  To  begin  let  R  denote  some  open 
connected  subset  of  sP.  The  set  R,  which  is  usually  determined  from  the 
physics  of  (1.1),  should  be  regarded  as  the  set  of  physically  admissible 
states.  We  first  require  that 

(2.1)  f(u)  f  C3( R)  , 

and  we  assume  that  the  eigenvalues  of  Df(u)  are  real,  distinct  and  are 
arranged  in  increasing  order: 

(2.2)  A  (U)  <  A  (u)  <  ...  <  A  (u)  . 

i  2  n 

In  addition  to  this,  we  require  that  for  each  A^(u)  and  ®H  u  e  R  either 

(GNL)  7uAk(u>  ’  rk(u)  ^  0  ' 

or 

(LD)  Vk(u)  *  Vu>  E  0  ' 

is  satisfied,  where  rk(u)  Is  the  right  eigenvector  of  Df(u)  corresponding 


to  X^(u);  see  [8].  Fields  that  satisfy  (GNL)  are  often  referred  to  as 
"genuinely  nonlinear"  fields,  and  fields  that  satisfy  (LD)  are  referred  to  as 
"linearly  degenerate". 

For  u_  (  R  it  is  well  known  that  (2.2)  implies  that  the  Hugoniot  locus 
is  locally  the  union  of  n  connected  one-dimensional  manifolds.  Following 
Mock  [11]  we  denote  the  Hugoniot  locus  about  u_  as 

n 

(2.3)  T(u  )  ■  U  r (u_)  , 

k-1 

where  the  index  k  is  chosen  such  that  the  tangent  vector  of  rk(u_)  at 
u_  is  rk(u_)«  We  require  that  the  local  structure  of  T(u_)  is  valid  in 
the  large  and  that  R  is  chosen  so  that  for  any  u_  f  R  we  have 

(2.4)  r<uj  C  R  . 

Next  we  impose  a  global  "entropy  condition"  on  those  fields  which  satisfy 
(GNL),  see  [8].  Specifically,  if  u+  c  rk(u_)  where  A^fu)  satisfies  (GNL), 
and  when  *k(u_)  >  *k(u+),  we  require  that 

(2.5)  Ak(u-}  >  s<u_'u+)  >  ' 

\+1<u+)  >  s(u_,u+)  >  # 

where  s(u_,u+)  is  the  "shock  speed"  determined  by  the  Rankine-Hugoniot 
condition.  (2.5)  is  of  course  the  celebrated  Lax  condition  E. 

Remark  2.1;  For  those  Xk(u)'s  that  satisfy  (GNL)  one  easily  deduces  that 
(2.5)  partitions  ^(u  )  into  two  disjoint  branches  separated  by  u_.  On  one 
branch  (2.5)  is  satisfied.  (Throughout  we  shall  denote  this  branch  by 
Ck(u_).)  On  the  other  branch  (2.5)  is  satisfied  with  the  roles  of  u_  and 
u+  reversed. 

Finally,  we  assume  that  f(u)  admits  a  convex  entropy  function.  A 
convex  entropy  function,  say  E(u),  is  a  scalar-valued  function  which 


satisfies: 


(2.6a)  E(u)  (  C3(R)  , 

(2.6b)  D^E(u)  is  positive  definite  on  R, 

5 

(2.6c)  D  E(u)Df (u)  is  symmetric. 

Given  that  n  ■  2  and  (2.2)  is  satisfied,  such  a  function  often  exists  [9]. 
Frequently  this  is  not  always  the  case  when  n  >  3.  Since,  as  will  be  seen 
below,  the  existence  of  E(u)  is  crucial  for  our  analysis,  we  should  mention 
that  numerous  physical  problems  do  indeed  admit  a  convex  entropy  function  [4] 
We  now  state  two  lemmas  without  proof •  Their  proofs  can  easily  be 
obtained  by  adapting  Lemma  5  and  Theorem  3  from  Mock  [11]. 

Lemma  2.1.  Suppose  that  all  of  the  assumptions  above  (with  the  possible 
exception  of  (2.6))  are  satisfied.  Let  u+  e  T^Cu^),  u+  /  u_,  where  A^tu) 
satisfies  (GNL) .  Then  the  function 

H(u_,u+ju)  =  -s(u_,u+)(u  -  u_ )  +  (f (u)  -  f(u_))  , 
has  no  zeros  in  I?  other  than  u_  and  u+.  Moreover,  given  u_  and  u+ 
and  any  open  set  0  containing  u_  and  u+,  there  is  a  positive  constant 
c(u_,u+,(?)  such  that  for  any  u  f  R \0 

| H(u_,u+;u) |  >  c(u_,u+,0)  . 

Lemma  2.2.  Given  that  E(u)  satisfies  (2.6),  there  exists  a  smooth  scalar¬ 
valued  function  F(u)  such  that 

(Df (u))TVE(u)  *  Vf(u)  . 

Moreover,  if  u+  u_)  lies  in  a  genuinely  nonlinear  field  on  which  the 
entropy  condition  (2.5)  is  satisfied  (i.e.  u+  c  T^(u_)),  then 
s(u.,u+) (E(u+)  -  E(u_) )  -  (F(u+)  -  F(u_ ) )  >  0  . 

On  the  other  hand,  if  u+  lies  in  a  linearly  degenerate  field  the  inequality 
above  becomes  an  equality. 


We  conclude  this  section  by  showing  that  the  flux  term  f(u)  from  the 


Navier-Stokes  equations  satisfies  all  of  the  assumptions  so  far  imposed. 
Recall  from  Section  1  that  the  flux  terms  of  the  Navier-Stokes  equations  are 
given  by: 


i 

I  f  (u) 


m 

m2/p  +  P 
m(e  +  P)/P 


where  P  is  the  density  of  the  fluid,  m  is  the  momentum  and  e  is  the 
total  energy  (kinetic  plus  internal  energy).  The  pressure  term  above  is  given 
by: 

P  -  (Y  -  1)(e  -  -i  m2/p)  , 

for  an  ideal,  calorically  perfect  gas,  where  Y  >  1  is  a  thermodynamic 
constant.  The  eigenvalues  of  Df(u)  are: 

(2.7)  1  ^  ( u )  =  m/P  -  C,  l2(u)  *  m/P#  *3(u)  *  “v/P  +  c  » 

where  C,  the  sound  speed,  is  given  by: 

C2  =  YP/P  . 

Furthermore,  A.j(u)  and  A3(u)  satisfies  (GNL),  while  X2(u)  satisfies 
(LD).  The  Hugoniot  locus  for  this  problem  is  readily  computed  to  be  [12]: 


(2.8) 


r1(u_) 


P+  *  P  0 


-1 


P+  =  P_(6  +  6)/(1  +  00) 


2/? 


( m/P ) +  =  ( m/p ) _  +  C_  _  1 


1  -  0 
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P+  =  P. 


r2(u_)  =  ^  p+  =  p_s 


(  m  /d  \  s  (  m  /n  \ 


>v>> 
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p+  -  P_e 


p+  -  p_(i  +  se)/(S  +  e) 


(m/P )+  ”  (m/p)_  +  C_  -  _  1 


8  -  1 
(1  +  60) 


1/2  ' 


"Y  +  1  y  -  1 

where  6  *  — -  ^  and  t  *  — — —  are  constants  and  the  parameter  9  lies  in 
0  >  0.  Moreover,  in  the  1  and  3  fields  above  the  entropy  condition  (2.5)  is 
satisfied  when  0  <  9  <  1;  see  Staoller's  book  [12]  for  a  thorough  development 
of  these  facts. 

We  take  for  the  set  of  physically  admissible  states  the  convex  set 
R  =  {(p,m,e)T  :  p  >  0,  P  >  0}  , 

and  one  easily  finds  that  T(u_)  C  R  for  any  u_  e  R.  It  is  furthermore  well 
known  that  the  ideal,  calorically  perfect  Navier-Stokes  equations  admit  a 


•.  ■.  i 


»v* 


.v.wv.' 

mi 

As  • 

AAA 


convex  entropy  function 


E(u)  -  -P  log  (P/( (Y-1 )P 1 ) )  , 


see  [4].  Clearly  E(u)  (  C3(R)  and  a  lengthy  calculation  will  verify  that 
E(u)  above  also  satisfies  (2.6b)  and  (2.6c)  [4,  1]. 

We  therefore  have  that  the  flux  terms  of  the  Navier-Stokes  equations  fit 
completely  into  the  framework  we  have  so  far  developed.  We  should  mention 
however  that  not  all  of  the  hypotheses  above  are  sacred  (other  than  (2.6)  that 
is).  We  have  chosen  the  route  of  simplicity  over  generality  here  since  such 
an  interesting  physical  example  is  so  easily  accommodated. 


§3.  LARGE-DATA  TRAVELLING -WAVE  SOLUTIONS  TO  POSITIVE  DEFINITE  SYSTEMS 

In  this  section  we  state  sufficient  hypotheses  and  outline  the  technique 
we  use  to  establish  the  existence  of  large-data  travelling-wave  solutions  to  a 
particular  version  of  (1.1).  The  specific  simplifying  assumption  we  make 


here,  in  addition  to  the  assumptions  of  Section  2,  is  that  the  "effective" 


V-V-Yj 

rm*-  ■  l 


diffusion  matrix,  D^EA,  has  positive  definite  symmetric  part.  Although  the 
main  result  of  this  section  is  a  generalization  of  previously  known  results 
(5,  2,  3,  11],  the  positivity  assumption  is  unreasonable  when  considering 
problems  such  as  the  compressible  Navier-Stokes  equations  or  other  physical 
examples  with  singular  viscosity  matrices.  Nevertheless  our  approach 
motivates  a  systematic  procedure  for  constructing  reasonable  artificial 
viscosity  terms  for  which  uniform  estimates  are  obtained.  With  these 
estimates  the  vanishing  artificial  viscosity  method  can  be  applied. 

Remark  3.1.  The  condition  that  D  EA  have  positive  definite  symmetric  part 
is  equivalent  to  the  condition  that  D^EA  have  positive  definite  symmetric 
part. 

We  use  the  method  of  continuation  to  obtain  a  large-data  travelling-wave 
solution.  A  homotopy  invariant  is  constructed  which  if  equal  to  unity 
guarantees  the  existence  of  a  travelling-wave.  TSie  problem  is  then  smoothly 
deformed  to  a  small-data  (weak-wave)  problem  for  which  it  is  routinely  seen 
(via  an  argument  similar  to  Foy  [5])  that  the  invariant  is  in  fact  unity. 
Before  stating  the  main  result  of  this  section  (Theorem  1)  we  outline  our 
method  of  proof. 

Let  Tk(u_)  denote  the  branch  of  the  Hugoniot  locus  on  which  the  entropy 
condition  (2.5)  is  satisfied,  (of  course  we  assume  that  X^(u)  satisfies 
( GNL ) ) .  We  suppose  that  u+  lies  on  this  branch  and  that  u+(a)  is  the  arc- 
length  parameterization  along  T^(u_)  with  u+(0)  =  u_  and  u+(L)  =  u+«  Now 
consider  the  variable  wave  strength  version  of  (1.1)  with  its  associated 
travelling-wave  equation: 


(3.1) 


A(w)  “  H(u_,u+(o) »w)  , 


lim  w(t)  ■  u_,  lim  w(t)  *  u+(o)  , 
t  +— 00  t++* 


where  H(u_ ,u+ (o) /u)  is  given  by; 

-s(u_,u+(a) ) (u  -  u_)  +  f(u)  -  f(u_)  . 

The  scalar  function 

(3.2)  A(u)  =  VE(u)  •  H(u_ ,u+ (c ) ;u) 

+  s(u_,u+(a))(E(u)  -  E ( u_ ) )  -  (F(u)  -  F(u_))  , 
will  serve  as  a  Lyapunov  function.  The  fact  that  A(u)  is  indeed  a  Lyapunov 
function  is  evident  by  differentiating  (3.2)  along  trajectories  of  (3.1), 
giving 


d  .  2  -1 

— —  A(w)  *  D  E( w)H(u  ,u  (C)w)  •  A  (w)H(u  ,u  (c);w)  . 
at  -  +  —  + 


(Recall  from  Lemma  2.2  that  DfTVE  =  VF.)  Since  by  assumption  D2ea  has 


positive  definite  symmetric  part.  Lemma  2.1  implies  that  A(w(t))  >  0  for 
all  w(t)  e  R  except  for  u_  and  u+(a).  Next  we  reparameterize  t  in 
(3.1)  so  that 

(3.3) 

Away  from  u_  and  u+(o)  this  is  equivalent  to  multiplying  the  left  hand 
side  of  (3.1)  by  the  positive  function 


A(w(t) )  =  1  . 


li(w)  «  D2E(w)H(u_,u+(a);w)  •  A  1  (w)H(u_,u+(a )  ;w) 


We  wish  to  show  that  the  unstable  manifold  of  (3.1)  near  u_  can  be 
connected  to  the  stable  manifold  near  u+.  With  this  in  mind,  consider  the 
region  in  state-space  given  by; 

(3.4)  M.=  LJt)  Oil.  , 

where 

L_(e)  *  {u  f  R  ;  A(u)  =  e),  (e  >  0)  , 

U  =  unstable  manifold  of  (3.1)  near  u_  . 


m 


,  * 


pp 


rm 


We  now  gives 


Claim  1.  For  all  0  <  a  <  L,  M_  is  homeomorphic  to  the  sphere  Sn"’*t  in  a 
small  neighborhood  of  u_  provided  that  e  >  0  is  chosen  sufficiently 
small.  (Note  that  throughout  this  section  k  represents  the  index  of  the 
entropy  condition  satisfying  branch  of  the  Hugoniot  curve  on  which  u+(a) 
lies . ) 

Before  proving  the  claim  we  state  a  lemma  from  Mock  [11]. 

Lemma  3. 1.  Let  B  and  C  be  matrices  and  suppose  that  B  has  positive 
definite  real  part,  Re(B)  =  1/2  (B  +  B*),  and  that  C  is  nonsingular  and 
Hermitian.  Then  the  form  q(z)  =  z*Cz  is  positive  definite  on  the 
generalized  eigenspace  of  B-1C  corresponding  to  those  eigenvalues  with 
positive  real  parts. 

Proof  of  Claim  1 :  The  proof  follows  from  the  assumptions  of  the  previous 

section  and  Lemma  3.1.  For  e  >  0  sufficiently  small,  L_(e)  is  a  manifold 

of  codimension  1.  This  is  clear  since  the  only  critical  values  of  A(u)  are 

A(u_)  *  0  and  A(u+(a))  >  0  (see  Lemma  2.2)  and  e  can  be  taken  between 

these  two  values.  The  unstable  manifold  of  (3.1)  is  given  locally  by  the 

generalized  eigenvectors  of  A  ^  (u_  )DH(u_  ,u+(cj  )  >u_ )  which  correspond  to  those 

eigenvalues  with  positive  real  parts.  Combining  the  facts  that  D*E(u_) 

symmetrizes  DH(u_ ,u+ ( a ) ; u_ )  and  D  EA  has  positive  definite  symmetric  part, 

it  is  easily  shown  that  the  number  of  eigenvalues  of  A  \u  )DH(u  ,u+(a);u  ) 

with  positive  real  parts  is  equal  to  the  number  of  positive  eigenvalues  of 

DH(u_,u+(a) ;u_) .  Counting  this  number  (that  is,  using  the  entropy  condition 

(2.5))  we  conclude  that  dim(U_)  =  n  -  k  +  1,  and  the  count  is  independent 

2 

of  0  <  a  <  L.  Finally,  using  Lemma  3.1  with  C  =  D  E(u_ )DH(u_ ,u+ (a ) ;u_ ) 

2 

and  B  =  D  E(u_)A(u_),  we  have  that 


‘ ^ V* V^T A^AT V, V V  V.  W  CTC«TT*g muwwswv  >"j  wjFjuvvMfij  ,r: 


'\F  vrwnr; 


is  positive  definite  when  acting  on  the  tangent  space  of  U_  near  u_. 

Thus,  M_  is  locally  homeomorphic  to  the  sphere  Sn_}t. 

What  follows  below  are  the  key  points  of  this  section.  Let  <j>t(x) 
represent  the  flow  of  the  system 

(3.5)  y (w)A(w)  ■  H(u_,u+(o);w)  , 

w(0)  -  x  , 

and  notice  that  <p.  (when  defined)  maps  M_  into  the  codimension  1  manifold: 
*0 


(3.6) 


L+(e)  «  {u  f  R  :  A(u)  ■  A(u+(o))  -  e}  , 


when  tg  =  A(u+(o))  -  2e  and  where  again  we  have  taken  e  >  0  small.  Next 
define 

(3.7)  M+  -  L+(e )  H  S+  , 

where 

S+  -  stable  manifold  of  (3.1)  near  u+(c)  . 

We  now  give  another  simple  claim: 

k—  1 

Claim  2.  For  all  0  <  cr  <  l,  M+  is  homeomorphic  to  the  sphere  S  in  a 

neighborhood  of  u+(o)  where  as  before  e  >  0  is  taken  sufficiently  small. 

The  homotopy  invariant  alluded  to  above  is  the  modulo-2  intersection 

number  of  <j>  (M  )  and  M.  defined  inside  the  ambient  space  L  (e);  see 
0 

[7].  As  per  Claim  1  and  Claim  2,  <p  (M  )  and  M+  have  complimentary 

0 

dimension,  that  is  dim(<p  (M  ))  +  dim(M+)  =  dim(L  (e))  ,  provided  that 

0 

<!>t  :  M_  ♦  L+(e)  is  a  smooth  map.  Unfortunately,  for  general  large-data 

C0 

problems,  this  need  not  be  the  case.  The  next  lemma  addresses  this  matter  in 
the  small. 

Lemma  3.2.  The  flow  <pt(M_),  0  <  t  <  tg,  remains  in  a  compact  subset 


of  R  for  all  0  <  o  <  an,  provided  that  o_ 


is  taken  sufficiently  small. 


Proof  of  Leona  3.2*  observe  that 

A(u+(o))  -  s<u_,u+(a) ) (E(u+(a) )  -  E(u_))  -  (F(u+(o))  -  F(u_))  . 

For  weak  and  entropy  satisfying  shocks  (i.e.  u+(o)  e  r^tu^)  with  a  small) 

Lax  has  shown  that 

A(u+(o) )  »  O(03)  >  0  . 

Computing  A(w(t(s))),  where  s  is  the  <.rc  length  parameter isat ion  along 
a  trajectory  of  (3.1),  we  obtain 

,  D2E(w)H(u  ,u(0)jw)  •  A  \w)H(u  ,u  (o)»w) 

—  A(w)  -  - — : -  * 

|A  (w)H(u_,U+(0) jw) | 

and  it  is  easily  seen  that  for  any  w  e  B(u^,6)  with  6  >  0  small,  we  have 

(3.8)  A(w)  >  const(1  -  <5)]h(u  ,u.(a)»w)|  . 

qg  "  t 

The  positive  constant  above  does  not  depend  on  6.  Expanding  Hock 'a  proof  of 
Lemma  2.1  one  determines  that 

(3.9)  |h(u_,u+(c) ?w) |  >  const  6c  , 

for  w  e  R  but  outside  B(u_,6)  U  B(u+(c),5),  again  for  6  >  0  and  o  >  0 
sufficiently  small.  We  next  show  that  w(t)  cannot  leave 
B(u_,p)  U  B(u+(a),p)  for  an  appropriate  choice  of  p. 

For  any  0  <  t  <  t  <  tg,  we  have  from  (3.3)  that 

t  a 

(3.10)  A(u+( 0) )  >  /  -jg  A(w)dt  , 

and  we  choose  t  to  be  the  first  time  that  w(t)  B(u_,5)  U  B(u+(c),6).  To 
reach  a  contradiction  suppose  that  w(t)  leaves 
B(u_,6  +  a)  U  b(u+(o),6  +  o).  Reparameterizing  the  right  hand  side  of  (3.10) 
in  terms  of  arc  length  and  using  (3.8)  and  (3.9)  we  must  have  that 

A(u+(c))  >  const(1  -  6  -  o)6o  •  o  • 


However,  recalling  that  A(u+(o))  =  O(a^),  the  inequality  above  would  imply 
that 

0(a)  >  const(1  -  6  -  a) 6  . 

Setting  6  *  ko  with  tc  <  1/(2a),  we  reach  an  obvious  contradiction  for 
small  o  and  large  k.  This  completes  the  proof  since  we  could  take  a 
smaller  if  necessary  to  force  B(u_,(ie  +  1 ) a )  U  B(u+(c),(ic  +  1)a)  to  be 
contained  in  R. 

Remark  3.2.  With  the  special  assumption  that  R  *  fP  along  with  the 

assumptions  that  D2E(u)A”^(u)  has  uniformly  positive  definite  symmetric  part 

and  A_1(u)  is  uniformly  bounded,  one  can  show  that  the  conclusion  of  the 

previous  lemma  remains  valid  independent  of  the  size  of  a. 

The  proof  of  the  next  proposition  follows  the  proof  of  Theorem  1. 

Proposition  1 «  Suppose  that  the  preliminary  assumptions  of  Section  2  are 

satisfied.  Furthermore,  assume  that  D2EA  has  positive  definite  symmetric 

part  in  R,  u_  e  R  and  u+(0)  €  r^(u^),  the  entropy  condition  satisfying 

branch  of  r  (u  ).  Then  I_(qp  (M  ),M ),  the  modulo-2  intersection  number, 
K  -  2  tQ 

equals  one  provided  that  0  <  o  <  aQ  with  aQ  sufficiently  small. 

Having  additional  knowledge  concerning  (M_ )  allows  us  to  state: 
Theorem  1.  Suppose  that  the  assumptions  of  Proposition  1  are  satisfied  along 
with  the  key  hypothesis: 

(A)  qpt(M_)  for  0  <  t  <  tp  remains  in  a  compact  subset 
of  R  for  all  0  <  a  <  L. 

Then  (1.1)  admits  a  smooth  travelling-wave  solution  for  any  u+(a)  with 
0  <  a  <  L. 

Proof  of  Theorem  1:  By  hypothesis  A,  Claim  1  and  Claim  2,  together  with  the 

usual  theorems  from  ordinary  differential  equations,  we  have  that  q>.  (M  ) 

0 

and  M+  are  compact  submanifolds  of  L.  (e).  Moreover,  qp. 

0 


and  M+  are 


smoothly  dependent  on  0  <  o  <  L,  and  they  have  complementary  dimensions  with 


respect  to  L+(e).  Therefore,  the  modulo-2  intersection  number, 

I-l?.  (M  ),M  ),  is  well  defined  and  remains  constant  throughout  the 
0 

entire  o  deformation.  Proposition  1  establishes  the  fact  that 

X2(<Pt  (M_),M+)  =*  1  for  small  o_,  thus  I-(et  (M_),M+)  *  1  for  any 
0  0 
an  <  o  <  L.  So  we  conclude  that  there  exists  at  least  one  trajectory  of 

'3.1),  (with  a  =  L),  that  connects  the  unstable  manifold  near  u_  to  the 

stable  manifold  near  u+(L).  This  is  the  desired  result. 

Proof  of  Proposition  Is  First  observe  that  (3.5)  can  be  written  as: 

(3.11)  U(w)D2E(w)A(w)  ^  =  7A(w)  . 

at 

Applying  the  results  of  this  section  it  is  clear  that  we  can  deform 
2 

D  E(w)A(w)  to  the  identity  without  changing  the  modulo-2  intersection  number, 
provided  0  <  a  <  oQ  and  oQ  is  sufficiently  small.  Now  let  {r^}”^ 
represent  the  orthonormal  set  of  eigenvectors  to  the  symmetric  matrix 

(3.12)  D2E(uJDH(u_,u+(cj);uJ  *  D2A(u_)  , 

arranged  so  that  its  eigenvalues  are  increasing  with  l.  Consider  the 

n 

orthogonal  change  of  coordinates  w  -  u_  =  £  v.r.  and  the  scalar  function 

l  1  1  1 


,c,  ,  1r~  ,2  .  b  ,  .3 

T  J,  W  *  *  <V 


and  b  can  be  rewritten  as 

b  -  -  (B^)2  +  0{o)  , 

where 

(8k)2  -  l?Xk(u_)  *  rk|(rjD2E(uJrk)  . 

(Again,  recall  that  the  index  k  refers  to  the  fact  that  u+(c)  e  rk ( u_ ) . ) 
Defining  v'  *  (v^#...#Vjt_^f0,vk+1,...#v  )  ,  we  have  from  Taylor's  Theorem 
that 

1^-  ( A(w )  -  Ac (v) )  -  0(|vk||v'|  +  |v'|2  +  |v|3)  , 

°K 

and 

( A(w)  -  Ac ( v ) )  -  0( | v | 2 )  for  l  t  k  . 

Rotate  coordinates,  (w  -  u_  =  Rv),  and  consider  the  final  deformation: 

dv 

p(Ti»v)  ~  *  7  A(n»v)  , 
at  v 

where  A(n;v)  =  Ac(v)  +  n(A(w)  -  AC(v))%  Following  the  outline  previously 

laid  down  in  this  section  (Claim  1,  Claim  2  and  Lemma  3.2)  along  with  the 

estimates  above  it  is  not  difficult  to  show  that  (M  ),M  )  is  well 

defined  and  remains  constant  for  all  n  e  [0,1]  provided  that  o  >  0  is 

fixed  and  is  sufficiently  small.  Setting  n  *  0  one  easily  verifies  that 

<p  (M  )  intersects  M+  exactly  once  and  that  this  intersection  is  indeed 
0 

transversal.  Therefore,  the  modulo-2  intersection  number  is  1  for  n  *  0 
from  which  by  homotopy  invariance  we  conclude  the  same  for  n  =  1. 

Remark  3.3.  We  have  intentionally  excluded  the  case  when  u+  lies  in  a 
linear  degenerate  field.  In  this  case  (3.1)  can  have  a  smooth  solution  only 
for  trivial  data  u+(o)  =  u_.  This  is  seen  by  first  recalling  that 


and  then  multiplying  (3.1)  by  D  E(w)i  doing  so  we  obtain: 

0  <  D2E(w)A(w)  ^  -  |r-  A(w)  . 
dt  dt  dt 

Lenina  2.2  states  that  when  u+(a)  is  in  a  linear  degenerate  field 
A(u+(a))  =  A(u_)  =  0,  therefore  the  equation  above  allows  us  to  draw  only 
one  conclusion,  that  is  w(t)  =  const. 

We  end  this  section  by  giving  sufficient  conditions  under  which 
hypothesis  A  of  Theorem  1  is  satisfied.  We  state  these  conditions  here  as  a 
theorem  since  it  is  this  technique  we  apply  in  the  following  sections. 

Theorem  2.  Let  G  C  if  represent  the  set  of  states  that  can  be  reached  by 
any  smooth  path  x(t)  with  x(0)  e  M_,  ^  A(x(t))  >  0  and 

e  <  A(x(t))  <  A(u+)  -  e.  Suppose  first  that 

(a)  G  C  R  . 

Second,  assume  that  there  exists  a  smooth,  nondecreasing  function  g(r),  with 
lim  g(r)  =  +“,  which  satisfies 

jr  +  ao 

(b)  min  (5Tli(u)D2E(u)A(u)5]  1/2  >  g’(|u|)  , 

Ul«i 

for  all  u  e  &  Then  remains  in  a  compact  subset  of  R  for  any 

0  <  t  <  A(u  )  -  2c. 

+ 

Proof:  Clearly  by  assumption  (a)  we  have  that  remains  in  a  closed 

subset  of  R .  To  see  that  <pt ( M_ )  remains  bounded,  let  m_  e  M_  and  observe 


[g( |w(t) | )  -  g( |m_ | ) ] 2  =  [/  g ( J w < t ) J ) dr ] 2  <  t  /  ||-  g( |w(t ) | ) |  dx 


Computing  the  derivative  in  the  right  hand  side  above  and  using  assumption  (b) 


we  arrive  at 


t  .  T 
dw 


[g(  |w(t)  |  )  -  g (  | m_ |  )  ]  2  <  t  /  -T~  li(w)D  E( w)A(  w)  —  dx  . 


Using  (3.11)  we  find  that  the  integrand  above  is  equal  to  A(w),  and 
recall  from  (3.3)  that  this  quantity  is  normalized  to  1.  Therefore  we 
conclude  that 

g(jv(t)|)  <  g(|m_|)  +  t  <  g(|ra_|)  +  A(u+)  -  2e  , 
which  implies  the  desired  result. 

$4.  SINGULAR  VISCOSITY  MATRICES 

We  carry  over  all  of  the  hypotheses  of  Section  2  to  investigate  problem 

(1.1)  in  the  case  when  the  viscosity  matrix  A  is  singular.  We  furthermore 
assume  the  set  R  of  physically  admissible  states  is  convex,  the  convexity 
of  R  along  with  the  assumption  that  D2e(u)  is  positive  definite  on  R 
allows  us  to  introduce  the  globally  defined  change  of  coordinates: 

(4.1)  v  -  VE(u)  . 

In  the  new  coordinates,  the  travelling-wave  differential  equation  (3.1) 
becomes : 

(4.2)  B( v)  ~  -  H(u_,u+»u(v) )  , 

lim  v ( t )  -  v^.  =  VE(u  )  , 
t*f 

where  B(v)  ■  A(u(v) ) (D2E(u(v) ) )” 1 .  Throughout,  B(v)  is  assumed  to  have 
positive  semi-definite  symmetric  part.  In  the  coordinate  system  given  by 
(4.1),  the  right  hand  side  of  (4.2)  is  given  by: 

H(u_,u+ju( v) )  -  VyA(u(v))  , 

(recall  that  A(u)  is  defined  in  equation  (3.2)).  To  simplify  the  notation 
below  we  define  ^(v)  by: 

A(v)  =  A(u( v) )  . 

Remark  4.1.  We  choose  here  to  work  in  the  coordinate  system  (4.1)  because  it 
is  this  coordinate  system  we  use  in  the  application  of  the  next  section. 


He  now  make  some  further  simplifying  assumptions.  Suppose  that  the  null' 


space  of  the  matrix  B(v)  is  spanned  by  p  independent  constant  vectors. 
This  allows  us  to  make  a  simple  rotation  (which  below  we  take  to  be  the 
identify)  so  that  B(v)  may  be  partitioned  as: 


(4.4) 


B(v) 


0  0 
b(v)  b(v) 


where  b(v)  is  a  (n  -  p)  x  (n  -  p)  positive  matrix  and  b(v)  is  an 
(n  -  p)  x  p  matrix.  To  further  simplify  our  presentation  we  assume  that 
b(v)  5  0,  and  we  note  that  the  case  when  b(v)  $  0  requires  only  a  slight 
modification  of  our  arguments  below. 

To  show  that  (4.2)  admits  a  solution  we  modify  the  diffusion  matrix  (4.4) 
by  introducing  an  artificial  viscosity  term.  Specifically,  we  consider  the 
family  of  modified  problems: 


(4.5) 


c,  e. 

B  (v  ) 


VvMve) 


lim  ve(t)  *  v  , 

t+±« 


where 


(4.6) 


Be(v) 


ea(v)I  0 

P 


b(v) 


and  where  a(v)  >0  is  to  be  determined.  Below  we  show  with  reasonable 
hypothesis  that  the  modified  problem  (4.5)  has  a  solution  for  any  e  >  0,  anc 
we  show  that  the  family  {ve}£>Q  satisfies  uniform  estimates  that  allows  for 
the  passage  to  the  vanishing  artificial  viscosity  limit. 

Before  giving  the  main  results  of  this  section  we  introduce  some  further 
notation.  Partition  a  vector  v  e  vP  as  v  *  (v1,V2)T,  where  v^  e  tP 
and  v2  e  Fn”p.  For  a  smooth  function  A  :  Vp  *  R,  let  A  represent  a 


vector  function  with  elements 


■* —  A,  1  <  i  <  p,  and  let  A  represent  a 

Yi  2 

vector  function  with  elements  w— —  A,  p  ♦  1  <  i  <  n.  The  notation  A 

1  32  .  11 

shall  represent  a  (p  *  p)  matrix  with  elements  gv  A,  1  <  i»j  <  p, 

i  ^  32 

and  Ay  y  shall  represent  an  (n  -  p)  *  p  matrix  with  elements  A, 

p  +  1  <  i  <  n,  1  <  j  <  p.  Finally,  for  any  square  matrix  M(v)  e  u*”*1"  we 


define  the  scalar-valued  function  Q,.(v)  by 

M 

Q^(v)  *  min  CTM(v)5  • 

|C  1-1 

We  are  now  ready  to  address  the  question  of  the  existence  of  solutions  to 
the  modified  problem  (4.5).  The  particular  list  of  assumptions  given  below 
are  chosen  for  convenience  only.  In  the  next  section  we  show  they  are  easily 
applied  to  the  compressible  Navier- Stokes  equations. 

Lemma  4. 1.  As  throughout,  suppose  that  all  of  the  basic  preliminary 
assumptions  of  Section  2  are  satisfied.  In  particular,  assume  that  u_  e  R 
and  u*.  f  r~  ( u_ ) ,  where  rfc(u_)  is  a  branch  of  the  Hugoniot  locus  on  which 
the  entropy  condition  (2.5)  is  satisfied.  Moreover,  assume  the  following: 

(a)  G'  C  VE(R )  , 

where  here  G'  C  a*1  is  the  set  of  states  that  can  be  reached  by  any  smooth 
path  y(t)  with  y(0)  •  v_,  ~  A(y(t))  >  0  and  0  <  A(y(t))  <  A(u+). 

Second,  assume  that  Q^tv)  >  0  and  a(v)  of  the  modified  viscosity  matrix 
(4.6)  is  constructed  so  as  to  satisfy 

(b)  c,  <  a(v)Q  .(v)  <  c2  , 

b 

for  some  positive  constants  Cj  and  c2  any  every  v  e  VE(R).  Finally, 
assume  that  for  all  v  e  G'  we  have 

(c)  [Qb(v)Q  ^(v)]172  >  g'Uelv^2  +  |v2|2)1/2)  , 

b 

for  some  smooth  nondecreasing  function  g(r)  with  lim  g(r)  =  +».  Then  for 


any  e  >  0,  the  modified  travelling-wave  equation  (4.5)  has  a  solution,  and 

g 

denoting  this  solution  by  v  (t),  we  have  the  estimate: 

g<(e|vE(t)|2  +  |vE(t)|2)1/2>  <  g(max(e,1)|v_|)  +  A(u+)  . 

Proof:  Observe  that  for  any  £  €  R11  we  have  that 


li(u(v))£TBE(v)£  >  |h(u_,u+»u)|2[j  Q  +  Q  .1(v)(ft2|2] 

b 


a 

I  2 


*  [Qa<v)|51|Z  +  Qb(v)|£2|2] 


where  Ift^2  +  |ft2|2  =  *  and  where  we  use  the  notation  £  *  (/F  £1#£2)T. 
Using  condition  (b)  one  easily  shows  that  the  right  hand  side  above  dominates 


|H(u_,u+ju)  |  2min(-^-,  iJmintc^Q^^tvJQ  ( v) )  | C | 


n2 


and  condition  (c)  gives  us  that  this  dominates 

c(|H(u_,u+iu)|g'(|v|)|£|]2  , 

for  some  positive  constant  c.  (Note  that  above  we  have  assumed  that 

g'(|v|)  <  c1#  of  course  we  lose  nothing  by  doing  this.)  Recalling  that 

Lemma  2.1  gives  us  that  |h(u_,u+ju)|  >  c(u_,u+,0)  for  u  outside  any  open 

dvE 

set  containing  u_  and  u+,  we  set  £  *  ^  and  mimic  the  proof  of 
Theorem  2.  This  completes  the  proof  of  the  lemma. 

The  estimate  of  Lemma  4.1  establishes  that  |vE(t)|  is  bounded 
independent  of  e ;  (provided  that  condition  (c)  does  not  depend  on  e.)  With 
an  additional  hypothesis  we  can  show  that  |vE(t)|  remains  bounded 
independent  of  e  and  therefore  obtain  a  uniform  maximum-norm  estimate. 

Lemma  4.2.  In  addition  to  the  hypotheses  of  Lemma  4.1,  assume  the  following: 
Suppose  there  exists  a  compact  set  I!  C  jf1  of  the  form: 

=  (v  :  v1  (  D  {v  :  |v2 |  <  R}  O  G' 
where  R  is  taken  large  enough  so  that 

g(R)  >  g(max(e,1) |v_| )  +  A(u+)  , 


[l 


R  R 

C  BP  contains  (v.j)_  and  (v1)+<  and  is  such  that  for  any 

v*  €  G 1  with  v*  f  3fl^,  | v* |  <  R  we  have  that 

K  (v*)  *  n  R 
1 

is  of  one  sign,  (n  t  Jp  is  the  outward  unit  normal  to  fl  . )  Then  v  (t) 
remains  trapped  in  (2  for  all  e  >  0. 

Proof :  By  the  estimate  of  the  previous  lemma  we  have  that 
g((e|v^(t)|2  +  | v| (t) | 2 > 1/2 )  <  g(R)  and  this  implies  | v^ 1 1 ) |  <  R.  Using 

the  differential  equations  (4.5)  and  the  assumption  of  this  lemma  we  have  that 

R  £ 

fi.j  is  either  forward  or  backward  invariant  with  respect  to  v^t).  Given 

£  R 

e  >  0  suppose  that  v^(t)  t  for  some  t.  If  this  were  the  case  then  we 

could  not  have  lim  ve(t)  ■  v  .  But  this  violates  the  result  of  Lemma  4.1 
t+±“  * 

and  therefore  establishes  the  result  of  the  present  lemma. 

The  next  lemma  implies  the  result  of  Lemma  4.2  as  well  as  giving  a 

dv£ 

uniform  maximum-norm  estimate  for  ^“* 

Lemma  4.3.  In  addition  to  the  hypotheses  of  Lemma  4.1,  suppose  we  could  find 

a  smooth  function  h  s  rP“P  ♦  tP  such  that  the  surface 

»  {(h(v2),v2>  :  v2  e  Rn"p}  O  G'  , 

contains  v_  and  v+  and  is  such  that  ^  (v)  ■  0,  (i.e. 

1 

Hj (u_,u+ ;u(v) )  ■  0),  for  every  v  (  Furthermore,  suppose  that  for  every 

v  e  E^,  the  (p  x  p)  matrix 

(<-s(u_,u+)I  +  Df(u(v)))(D2E(u(v)))*1]i>j  1  <  i,  j  <  p  , 

(■A  (v)) 

V  V 
11 

is  either  positive  or  negative  definite.  Then  ve(t)  satisfies  the  uniform 


estimates: 


«  K1  , 


-•  -»  -  *  -  •-  -  1-mA  - - 1 - -A-  A  -  lLA  %  -  «•-  A.  »•  V  V  S  •  J*  -  >  >  k  •- 


>dt  ■-  '  “2  ' 

for  some  positive  constants  K^  and  Kj. 

Proof:  He  first  prove  (a).  The  idea  is  to  show  that  ve(t)  stays  in  a  small 

neighborhood  of  E1  for  0  <  e  <  eQ  with  eQ  sufficiently  small.  Set 

eJ  «  {v  :  dlv,^)  <  6}  , 

with  6  small  but  fixed  and  define 

E*e  -  {v  :  | A  < v) |  <  a(v)Ke>  . 

1 

Defining  R  as  in  the  previous  lemma  choose  Ke  sufficiently  small  so  that  a 
connected  branch  of 


Ke 


n  {v  :  I v_ I  <  R}  n  G' 


is  contained  in  E  .  Next,  let  n  represent  the  outward  unit  normal  of 


E*0,  and  compute  that  along  SE^  O  {v  s  |v2|  <  R}  Pi  G  • ,  n  is  given  by 

| A  (A  ,0)T|“1»(A  A  ,A  A  )T  . 

V1  11  1  v2  1  1 


Ke 


Dotting  (4.5)  with  n  gives 


dv 

dt 


•  n  -  |A  (A  f0)T|"1|A  | 

1  trtr  it  1  1  If  1 


w  v. 


l"Av  I 

r  T-  1  -1- 

[C  A  c - —  +  (b  A  ) 

<-  v  v  ca  v 


V1V1  ^ 


”1’  )TA  s]  , 
v2  vzVl  J 


where  £  *  A  /|A  |.  Since  along  3E  we  have  |A  |  =  aKe,  and  since  the 

V1  V1  _  V1 
hypotheses  of  the  lemma  gives  us  that  A  is  either  positive  or  negative 

11 

definite  near  E^,  we  have  that  the  bracketed  term  above  can  be  made  to  have 
one  sign.  Itiis  is  accomplished  for  all  v  e  3E2E  n  {v  :  |v2|  <  R}  n  G’  by 
choosing  K  sufficiently  large  when  e  is  sufficiently  small.  Following  the 
reasoning  of  the  proof  of  Lemma  4.2  completes  the  proof  of  (a). 

To  prove  (b)  observe  that 


I  A.,  tv)  I 


Since  v£(t)  e  Z*e  n  {v  s  |  |  <  R}  HG*  we  have  |~  v^(t)  j  <  K  for  all 
sufficiently  small  e  and  all  t.  Combining  this  with  the  trivial  inequality 
(t) |  K  jb  1(v)Ay  (v)|,  completes  the  proof  of  (b). 

We  are  now  ready  to  state  and  prove  the  main  result  of  this  section. 
Theorem  3.  Suppose  that  the  assumptions  of  Lemma  4.1  and  Lemma  4.3  are 
satisfied.  Then  the  partial  differential  equation  (1.1)  with  singular 
viscosity  has  a  smooth  travelling-wave  solution  u(x  -  st)  which  is  the  limit 
of  a  sequence  of  artificial  viscosity  approximations. 

Proof:  From  Lemma  4.1  we  have  a  family  of  smooth  solutions  to  (4.5C)  which  we 
denote  by  {v  }e>Q*  Normalize  these  by  translating  t  so  that 

A(ve(0 } )  -  -1  A(u+)  . 

Moreover,  from  Lemma  4.3  and  the  Arzela-Ascoli  Theorem  we  have  a  continuous 

ek 

function  v(t)  and  a  subsequence  v  (t)  such  that 

ek 

v(t)  -  lirn  V  (t)  , 


the  convergence  being  uniform  on  compact  t  intervals.  From  the  proof  of 
Lemma  4.3  we  also  have  that 


Hi (u_,u+ju( v(t) ) )  *  lim  A  (v  (t) )  »  0 


e.+O  V1 
k 


and  the  usual  bootstrap  arguments  from  the  theory  of  ordinary  differential 
equations  gives  us  that  v2(t)  is  a  smooth  solution  of 


b(v)  -  H2(u_,u+?u(v)) 


Since  by  assumption  A  is  nonsingular  for  v  e  Z.,  we  have  from  the 

11 

implicit  function  theorem  that  v^(t)  is  smooth  as  well.  Therefore,  v(t) 


-23- 


is  a  smooth  solution  of 


B°(v)  ~  =  H(u_,u+;u( v) )  . 


What  remains  to  be  shown  is  that 
compute  that 


lim  v(t) 
t+±® 


v^.  With  this  in  mind 


(4.7) 


4-  *<v(t))  =  l  ( v( t) ) 
dt  v2 


b-1(v(t))K  (v(t))  . 

v2 


Therefore,  lt(v(t) )  is  a  nondecreasing  function,  and  because  it  is  the  limit 

_  ek 

of  A(v  (t)),  which  is  bounded  between  0  and  A(u+),  we  find  that 

0  <  J(v(t) )  <  A(u+)  . 

Bounded  monotonic  sequences  have  limits,  therefore 

lim  A( v(t) )  =  A  , 
t>±“ 

exist.  Now  define  v^n^(t)  by 

v(n)(t)  =  v(t  +  n)  , 

and  again  appealing  to  the  Arzela-Ascoli  Theorem,  we  have  a  continuous 

function  v*(t)  and  a  subsequence  nk  ♦  ®,  such  that 

(nk) 

v*(t)  =  lim  v  (t)  , 

V 

which  also  gives  us  that 

A  =  lim  A(v(t  +  n,  ) )  =  A(v*(t))  • 

+  k 

n  +oo 

k 

Integrating  (4.7)  from  t  *  n^  to  t  =  1  +  n^,  we  have  after  changing 
variables 

_  (v  .  (V 

A ( v  (1))  -  A( v  (0)) 


1  (nk)  (nk>  -  (V 

=  /  A  (V  (t))  •  b  (V  *  ( t ) ) A  (v  (t))dt  . 


Letting  +  •  implies  that 

0  -A  (v*(t))  •  b"1(v*(t))A  (v*(t))  , 

2  2 

and  therefore 

X  (v*(t))  ■  H2(u_,u+»u(v*(t) ) )  -  0  . 


Combining  this  with  the  fact  that  H.|(u-#u+;u(v*(t) ) )  “0,  we  have  that 

u(v*(t))  is  one  of  the  two  critical  points  of  H(u.,u+;u).  However  because  of 

the  normalization,  that  is  A(v(0))  >  —■  A(u+),  we  must  have  that 

v* (t )  =  v+.  The  same  argument  can  be  applied  equally  well  to  any  subsequence 

of  v^n^(t)  ■  v(t  +  n)  which  shows  that  lim  v(t)  *  v+.  Similarly, 

t+*» 

lim  v(t)  *  v.  and  so  the  proof  of  the  theorem  is  complete • 
t>— <® 


§5.  APPLICATION  TO  THE  NAVIER-STOKES  EQUATIONS 

In  this  section  we  apply  Theorem  3  to  the  compressible  Navier-Stokes 
equations •  From  our  analysis  we  conclude  that  these  equations  admit  a  smooth 
travelling-wave  which  is  the  limit  of  certain  artificial  viscosity 
approximations.  Throughout  this  section  we  shall  assume  that 


u 


{(p,m,e) 


T 


p  >  0,  P  >  0}  , 


and  we  assume  that  u+  e  r^(u_)  where  k  «  1  or  k  *  3  (the  genuinely 
nonlinear  fields);  see  equation  (2.8).  Moreover,  because  the  Navier-Stokes 
equations  are  invariant  under  a  Galilean  change  of  coordinates,  we  lose  no 
generality  by  taking  the  shock  speed  s  equal  to  zero. 


We  begin  by  explicitly  transforming  the  travelling-wave  equations  for 
Navier-Stokes  into  the  coordinate  system  given  by  the  gradient  of  its  entropy 


function 


E(u)  «  -  P  log(p/((Y-1>PY))  . 

Computing  some  derivatives  we  find  that 

vt  =  0  =  -log(2pe  -  m2)  -  2pe/(2pe  -  m2)  +  (y  +  1)(log(p)  + 

v2  =  |^  =  2pm/(2pe  -  m2)  (-(Y  -  1)  ^)  # 

v3  =  ||  =  “2p2/ ( 2pe  -  m2)  (=  -(Y  -  D  f)  • 


Inverting  these  equations  we  obtain 

-1 

p(v)  =  ( -v 3 ) ^ " 1  exp [—  j (v1  +  v2/(-2v3>  -  Y)]  • 
m  =  p(v)v2/(-v3)  , 
e  «  p(v)(1  +  v2/(-2v3) )/(-v3)  . 

Therefore 

f-j  *  m  «  p(v)v2/(-v3)  , 

2  2 

fy  -  —  +  p  =  P(V)(1  -  Y  +  V  /v )/v  , 

*■  p  2  3  3 

f3  =  p  (e  +  P)  “  P(v>v2(y/v3  "  \  (V2/V3)2^/V3  ' 


and  it  is  easily  verified  that 

f(u)  -  f(u_)  =  ^(v)  , 

where 


A(v)  =  (y  -  1)P(v)v2/(-v3)  -  f(u_)  •  v  +  const 
and  where  the  constant  is  chosen  so  that  A(v_)  *  0.  Moreover, 
lengthy  calculation  will  reveal  that 


1)  +  log  2 


a  rather 


B(v)  5  A(u(v) ) (D2E(u(v) ) )_1 

0  0  0 
0  -1/v3  v2/v23 

0  v2/Tj  -v2/vj  ♦  1/v 

V 

(The  viscosity  matrix  A(u)  for  Navier-Stokes  is  given  in  Section  1.)  Note 
that  in  the  v-coordinate  system  the  region  of  physically  admissible  states  is 
VE(R)  -  {v  :  v3  <  0}. 

The  result  of  this  section  is  obtained  by  establishing  the  hypotheses  of 
Theorem  3.  To  this  end  we  first  establish  hypothesis  (a)  of  Lemma  4.1.  That 
is  we  show  that  the  set  6'  is  contained  in  a  closed  subset  of  VE(R). 
(Recall  that  G*  is  the  set  of  all  states  that  can  be  reached  by  any  smooth 

path  y(t)  with  y(0)  -  v„,  ^  X(y(t))  >  0  and  which  satisfies 

0  <  I(y(t))  <  7(v+).)  We  accomplish  this  by  constructing  what  we  call  a 
"A-wall".  Specifically,  we  show  that  when  k  *  1,  m_  >  0  (The  case 
k  »  3,  m_  <  0  being  similar)  there  exists  a  closed  set  of  the  form 

111  »  {v  :  Vj  )  a  >  0,  Vj  <  B  <  0>  , 
with  v_  e  (0°,  v+  e  3(i I  such  that  for  any  v*  e  3W\{v+}  we  have 
J5(v*)  >  J(v+).  Having  this  wall  Ul  implies  that  G'  C  VE(R). 

Remark  5.1.  The  entropy  condition  (2.5)  implies  that  the  momentum  component 
of  u»,  (i.e.  m_)  for  the  Navier-Stokes  1  or  3  zero-speed-wave  can  never  be 

zero.  See  (2.7)  and  compare  it  with  the  entropy  condition  (2.5). 

To  see  that  U/  can  be  made  to  have  the  properties  described  above,  set 

a  •  (v2)+  , 

&  -  -4(y  -  1)/(f2(u_)/mJ2  . 


*  > 


vXw 


“•  ,v 


v\-> 

vVv 


and 


(From  the  remark  above  (f2(u_)/m_)^  is  positive  and  bounded.)  By  checking 

g 

the  sign  of  -r— —  X(v*)  for 
3 

v*  €  {v  s  H(v)  =0,  v2  -  (v2)+,  v3  <  0}  , 


it  is  clear  that  H(v*)  >  A(v+)  when  v*  ?  v+;  see  Figure  1.  Moreover,  it 

g 

is  easy  to  check  that  X(v*)  >  0  for 


/*  e  {v  s  X(v)  =  0,  v2  >  (v2)+  ,  v3 


4(y  -  1)  1 

2>  ' 

(f2(u_)/m_) 


so  we  have  A(v*)  >  A(v+)  here  as  well;  again  see  Figure  1.  Finally, 

H(v*)  is  positive  (resp.  negative)  if  v*  c  3W  with  v*  lying  above 

1  A 

(resp.  below)  the  surface  {v  :  ^  ■  X(v)  *  0,  v2  >  0,  v3  <  0},  so  we  conclude 

that  A(v*)  >  A(v+)  for  all  v*  e  3W\{v+}. 


Figure  t 


Next  we  establish  hypotheses  (b)  and  (c)  of  I^mma  4.1.  Define  the 


viscosity  submatrix 


and  construct  the  artificial  viscosity  ea(v)  (■  ^(v))  so  that  for  all 

v  t  {v  :  v3  <  0} 


c-  <  a(v)fi  .(v)  <  c2  , 
b 

where  c1  and  c2  are  some  positive  constants.  Note  that  for  a  symmetric 
matrix  S,  Qg  is  given  by  the  smallest  eigenvalue  of  S. 

Remark  5. 2.  This  particular  artificial  viscosity  has  the  form  (1.5)  when 
written  in  conserved  variables. 

An  elementary  calculation  will  show  that 


Qb(v,Q  -1(v) 


and  for  any  v  e  U  this  dominates 

2 

c 

(1  +  |v2|  +  |v3|)2 

for  some  positive  constant  c.  With  this  estimate  we  can  satisfy  hypothesis 
(c)  of  Lemma  4.1  by  choosing  g(r)  *  c  log( 1  +  r). 

To  conclude  this  section  we  show  that  the  Navier-Stokes  equations  satisfy 
the  assumptions  of  Lemma  4.3.  Solving  the  equation  fi(v)  «  0  for  v1  we 


find  that 


v ^  -  h(v2#v3) 


(-2v3) 


—v  y— 1 

+  Y  +  (Y  -  1)log((— 21  m  )  , 

2 


and  clearly  this  is  smooth  in  W.  Finally  for  v  *  (h(v2#V3),V2»v3),  that  is 

for  v  such  that  t— —  J(v)  -  0,  we  have  that 

<jv  . 


A(v) 


3v‘ 


(Y  “  D 


and  as  Remark  5.1  points  out,  m_  can  never  be  zero  for  a  zero  speed  1  or  3 
wave.  Therefore  we  have  established  all  of  the  hypotheses  of  Theorem  3  and 


«VV. 


hence  conclude  the  result  of  Theorem  3  for  the  Navier-Stokes  equations 
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